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Abstract 

We report on a calculation of the effects of supersymmetry on the cross-section 
for t — i production at LHC. A numerical study is carried out for the ten bench- 
marks of the Snowmass accord. It is found that the higher order effects involving 
supersymmetric particles in internal loops can be as high as 6%, both for the 
cross-section and the (parity even) helicity asymmetry, for one particular bench- 
mark. For other benchmarks smaller but nonetheless observable corrections are 
found. 
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1 Introduction 



For nearly 35 years, supersymmetry (SUSY) has been an attractive theory in particle 
physics. At the most theoretical level, it permits the construction of string theories 
which do not contain tachyonic states, and at the phenomenological level it offers an 
explanation for the naturalness of the hierarchy through its reduced ultraviolet diver- 
gences, as well as providing resolutions of several puzzles arising in standard models of 
cosmology. It also gives rise to a correction to the running of the couphngs, so that the 
strong, weak, and electromagnetic interactions can unify at some Grand Unified (GUT) 
scale. 

However, to date there has been no reliable evidence that this theory describes 
Nature, so that if SUSY is indeed realised in nature, it must be broken at a scale higher 
than that reached in accelerator experiments conducted up to now. If the theory is to 
be effective in providing a solution to the naturalness of the hierarchy problem, then 
the supersymmetry breaking scale cannot be much more than about 1 TeV. This is also 
the scale of SUSY breaking which leads to unification of couplings. This implies that, 
with the exception of some hidden corners of parameter space, SUSY is expected to be 
discovered at the forthcoming LHC. 

Clearly the most dramatic manifestation of SUSY would be the production and iden- 
tification of super symmetric partner particles such as the spin-| charginos or neutralinos, 
or evidence that at sufficiently high energies hadrons display behaviour consistent with 
the existence of squarks or gluinos in the sea. Nevertheless, the existence of supersym- 
metry will have indirect but measurable effects on the cross-sections for the production 
of Standard Model (SM) particles. The LHC is expected to achieve sufficient integrated 
luminosity such that it will be possible to determine these cross-sections with sufficient 
accuracy to be able to detect the effects of higher order corrections coming from loops 
of supersymmetric particles. The loops can give rise to a significant correction to the 
production cross-sections even below the threshold for the production of the supersym- 
metry particles themselves, so that hints that some new physics is imminent can be 
deduced before the threshold energies are actually reached. Above these thresholds, dif- 
ferential cross-sections with respect to suitably chosen variables can display structures 
which can be used to determine masses of some of the scalar particles of the MSSM. 

In this paper, we consider the influence of the minimal supersymmetric extension 
of the Standard Model (MSSM) to the production of top-antitop t — t. At LHC, after 
an integrated luminosity of 100 fb~^, one expects around 2 x 10^ such events, so that 
small corrections should be easily identifiable, not only at the level of the total produc- 
tion cross-section, but also for differential cross-sections with respect to the transverse 
momentum, pt of the t-quark, and also with respect to the invariant mass Mtt, of the 
t — t pair. Both of these are expected to be measured accurately at LHC. The latter 
variable is one for which differential cross-sections display a discernible structure as one 
crosses various thresholds. Apart from statistical errors, there will be larger systematic 
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errors arising form uncertainties in the incident beam flux and in the parton distribution 
functions (PDF's). Such systematic errors are cancelled in the ratio of helicity asym- 
metries to total cross-section, defined in Eqs. (14.21) and fl4.3l) . Although the helicities of 
individual t- or t-quarks cannot be determined on an event-by-event basis, their distri- 
butions can be inferred from the angular distributions of the decay products (see [1]). 
The corrections to these asymmetries can be as high as 5 %. 

Although the MSSM has very few extra parameters in the supersymmetry conserv- 
ing sector, the supersymmeric breaking terms can introduce 105 parameters in addition 
to the 19 parameters of the SM. Clearly, it is impossible to analyse the complete space 
of these parameters. Nevertheless, we have therefore organised our calculation so that 
software is available to calculate differential cross-sections with all possible helicity con- 
figurations and any given set of SUSY parameters. We have done this by working in 
terms of helicity matrix-elements and setting up "prototype graphs" which can then be 
included in the calculation with any combination of couplings and internal masses. This 
provides maximum flexibility for the extraction of total or differential cross-sections, he- 
licity asymmetries, etc. as well adaptability to other SUSY models. The library was also 
designed to scan efficiently over a large number of parameter sets. This was achieved by 
keeping, for one set of kinematic parameters, a table of all required Veltman-Passarino 
functions in memory, thus avoiding to re-calculate these functions for each new set of 
model parameters. This can reduce the computation time by up to 60% and makes 
scans over one or even two-dimensional parameter spaces feasible. 

The number of independent SUSY parameters is greatly reduced if one considers 
supersymmetry scenarios which are consistent with super-gravity or in which SUSY 
breaking is either gauge-mediated or mediated by the Weyl anomaly. At the Snowmass 
meeting of 2002, [5J, a collection of ten typical models was considered and sets of SUSY 
parameters for these models were generated. For convenience, we have investigated the 
effects of SUSY corrections to t — t production for these ten parameter sets. We find 
considerable variation in the magnitudes of the corrections from these different param- 
eter sets. Conversely, this means that accurate measurement of the t — i production 
cross-section can be used as a tool to help identify the correct set of SUSY parameters. 

At sufficiently high (partonic) energies, the SUSY corrections to t — i production are 
expected to be dominated by single and double logarithms of incoming parton energy 
divided by the SUSY breaking scale, Msusy- The determination of these logarithms is 
independent of the SUSY parameter set, with the exception of Msusy and the ratio, 
tan/3, of the vacuum expectation values of the two Higgs doublets, and the calculation 
is simplified by the fact that the mixing of various supersymmetric particles to form 
mass eigenstates has no effect on these logarithms. The logarithmic contributions have 
been calculated by Beccaria et. al. [5]. One may have expected that it would have 
been possible to express the entire SUSY correction in terms of these logarithms plus 
a constant off-set, which depended on the SUSY parameter set. We have compared 
our results with those of ref. [9] and although it is indeed the case that our results 
agree with these logarithms plus a constant off-set at sufficiently high partonic energies, 
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this approximation is found to be unsuitable at typical partonic energies which will be 
reached at LHC, and the entire calculation is required for a reliable prediction of the 
cross-sections at LHC. 

The structure of this paper is as follows: In section^we discuss the general method 
for the extraction of the above-mentioned helicity matrix-elements at the partonic level 
from a general Feynman graph. In section^we list all the prototype graphs and indicate 
which supersymmetric particles can contribute for each of the prototypes. In section^ 
we discuss the results after folding the partonic cross-sections with PDF's and show the 
results for the ten Snowmass benchmark points, ^ect^on [5] presents some conclusions. 



2 Helicity Matrix Elements 

Because we wish to be able to discuss the total and differential cross-sections for given 
helicities of the t- and t- quarks as well as the asymmetries, we find it convenient to 
work at the parton level for a generic process 

in terms of helicity matrix-elements ^Aa,A6,Ai,A2(-^7 ^)) functions of the partonic centre- 
of-mass energy, E and scattering angle 6 (also in the centre-of-mass) . 

In the case of quark-antiquark annihilation, for which we take the incoming quarks 
and antiquarks to be massless, the helicities of the incoming partons are anti-correlated 
(Aq = —Aft), and although it is possible for this initial helicity anti-correlation to be 
violated at the one-loop level, the interferences at 0{ag) or O^a^aw) always respect 
this anti-correlation. On the other hand, for the gluon fusion process there need be no 
correlation between gluon helicities. 

Once these helicity matrices have been determined, differential cross-sections and 
asymmetries can be computed by convolution with the corresponding parton distribu- 
tion functions (PDF), summed over helicities or not, as appropriate. 

A further potential advantage of the helicity matrix-element formalism, although not 
currently applicable at LHC, is the determination of initial beam polarisation asymme- 
tries, should it become possible in the future to polarise these beams. For interactions 
of a parity violating nature such as SUSY corrections to weak interactions, such asym- 
metries would be immensely useful in identifying the parameters of the supersymmtery 
model. 

The helicity amplitudes are obtained in two stages. In the first stage, a set of 
coefficient functions, a\"l^{E,6), of complete set of Dirac 7-matrices, is determined: 

A,MMiE,e) = Y.at:\^{E,e)u{p,,x,)r^-^v{p2,-X2) (2.1) 

{a} 



4 



where T^"^ are the matrices 



V 


= r 


J- A 


= i^i 







(2.2) 



with associated projection operators P't"^ 



— 




^A — 




JTj, — 


c 



5 Pi 



(2.3) 
, such that 

Tr(r^">P{,}) = 6\;l. 
Note that the coefficient functions are independent of the hehcities of the t- and t quarks 
For the basis vectors Cq ■ ■ ■ 63, where 

= ^(-(p + ^cos%^ + (p-Ecos%^ + 2pp^) 

^ 2E2psin/^'"^^'^^Pa 

4 = ^iPl-P2) 

{p being the magnitude of the three-momentum of the t-quark in the centre-of-mass 
frame), the hehcity matrix elements are given by0 

A = 2[E {—ay + iXitty) + ip (a^ + iXia\) — irrit [a^ — iXia^)] ^Xi-\2 

+ 2 [mta^ -pa^2 -iEAia?? -mtAia^] 5ai,A2 (2.4) 

Thus, for example, the non-zero coefficient functions for the quark- ant iquark anni- 
hilation process at the tree-level (Figure [1]) are: 

T 2 cos 6^ ^ , , 

«y = ^.-^^a^-aJt^t) 

4 = ^?'^'^a„,-aJt ® r) 

1 asin 6* ^ , , ,^ 

4 = 9s^Sxa,~xAr®T), (2.5) 



^ These matrix-elements are defined up to an overall phase, which may depend on the initial-state 
and final-state helicities. 
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Pa(Aa) 

Figure 1: Tree-level graph for t — t production from quark-antiquark annihilation. 




where (r ® r) indicates the colour factor for a single gluon exchange. Inserting these 
expressions for the coefficents into eq. fl2.4p generates the helicity matrix-element 



-9I (Aa + cos^) (5a„,-a6^Ai,-A2 (t 
-9s ^ ^a,,-A6<3ai,a2 (t O 1 



(2.6) 




.-Q-Q-O-O-Q-Q. 




Figure 2: Tree-level graphs for t — t production from gluon fusion, 
whereas for the gluon fusion process (Figure [2]), the non-zero coefficients are given 



by 



9! 



{E sin 9 (cos 9-2) 5\^^x^^ + 2p sin 9 cos 9} 



2E{E-pcos 
+ [r'^ ^ A, ^ A,, (tt + ^)] 



2E{E-pcos( 



,12 



-rV> sin ^ A,5a„,-a, + [r'^ r^ A, <^ A5, ^ ^ (vr + 9)] 
-t''t^ { ~E cos 9 (cos ^ - 2) 5a„,a, + 2p sin^ ^ } 



2E(E -pcos6 

r\ A, ^ A5, e ^ (tt + 



-rVA, (pcos^-F)5a., 



2E(E -pcos( 
+ [r" ^ r^ A, ^ A,, ^ (tt + 



2E{E-pcost 



-/r^mt cos^(5a„,a, + [r^ r^ A, ^ A,, ^ (vr + 9)] (2.7) 
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leading to a tree-level helicity matrix-element 



-r T 



2E(E -pcosi 

X \mt [2Ai {{E + p) cos' d - 2p ~ AE cos d) - 2\aE] Sx,m^x^,> 

[-2E{E + p) cos ^ sin ^ + AE^ sin 9] 5a„-a,5a„,a, 
-2\ipmt sin' 6'5ai,a2'^a„ -a, 

-£;psin6' (2cos6' + AiAa) 5ai -Aa^A^ -a^ 

+ [t'^ ^ t\ Xa ^Xb, e^{7r + 6)] (2.8) 

In the case of gluon fusion the contribution from any graph may be written in the 
form 

m(Pi, Xi)Tv{p2, -A2), 

where F is a sum of strings of 7-matrices with coefficients that are proportional to 
couplings, internal and external fermion masses and the Veltman-Passarino (VP) [2] 
functions arising from the loop integrals. These VP functions have arguments that 
depend on the internal and external masses as well as on the Mandelstam variables 
s, t, u. The coefficients a^"'^ are simply projected by 

= Tr (P^">r) . (2.9) 

For the quark-antiquark annihilation process there are two types of contributing 
graphs: 

The first type are graphs for which the fermion lines can be factorised into an initial 
quark line, F' and a final t-quark line, Tf. Again these are sums of string of 7-matrices 
with coefficients that are proportional to VP functions. In this case the coefficients a^""^ 
are projected by 



^A„,-A,Tr (va ■ ^^^^^) Tr (P^">F;) , (2.10) 



where v^^ is a vector in the plane normal to the incoming momenta Pa and pt, given by: 

1 



{Pb ■ PlPa + Pa ■ PlPb - Pa " PbPl + i^aS'^paPaPhPl} 



^ /?: 1-1 ■ n L-f'" fi-t'a ' fa fi-fb t^a t'uf i ' "-'a^i/pat 

y 2Ep SlYlU 

The other type of graph is one in which the fermion lines do not factorise into an 
initial fermion line and a final-fermion line, but rather into an upper line, F„, and a 
lower line, F^, in which the incoming quark and outgoing t-quarks are connected by the 
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Pa{K) ► 



Pl(Al) 



?'2(A2) 



Figure 3: Gluino, neutralino, or chargino exchange contribution to t — t production from 
quark-antiquark annihilation 

exchange gluino, neutralinos, and charginos,is the s-channel. It is graphs of this type 
that give rise to non-zero amplitudes in the case where the incoming helicities are equal, 
but such amplitudes do not interfere with the tree-level amplitudes. An example of such 
a graph is shown in Figure [31 

For such graphs the coefficients a^"^ are projected by 



As in the case of most of the contributing graphs, the expressions obtained from 
FigureOfor the corrections to the coefficients are too long and unwieldy to be reproduced 
here 



Figure 4: Graph involving the exchange of intermediate Majorana fermions 

Finally we note here that in several cases, the internal fermions exchanged in the 
s-channel in such graphs may be neutralinos or gluinos, which are Majorana fermions. 
In such cases supplementary graphs of the type shown in Figure H] need to be considered. 
Great care needs to be taken in handling such graphs. The standard expressions for 
the propagators of Majorana fermions in which a fermion propagates into a fermion or 
an ant-fermion propagates into and anti-fermion, are ambiguous up to a sign until the 
exact ordering of the fermions is determined. In order to ensure that this is effected in 
a consistent manner it is necessary to determine the fermion ordering of the term in the 
Wick contraction corresponding to the graph under consideration. 

A library has been constructed both in FORTRAN and C-I--I- in which each of the 
prototype graphs shown in the next section can be determined numerically, as a function 
of the incoming energy, scattering angle, helicities, couplings and internal masses. We 
have checked all prototype graphs by selecting different routings of the internal loop 
momenta and copmparing the results numerically. The numerical values of of the rel- 
evant VP functions are determined either using the FF library [3] or LOOPTOOLs [1] 
These libraries can be found at 

http:/ /hep. phys.soton.ac.uk/ /hepwww/staff/D.Ross/susyttbar/ 




(2.11) 
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3 Diagrams 



In this section we list the diagrams needed for the computation of SUSY contributions 
to polarised top-antitop production cross sections. To save space we only draw the 
different topologies. For further reference, each topology has a label set in typewriter 
font. An asterisk behind the label indicates that the crossed version of this diagram 
has to be included as well. A double asterisk indicates that only the crossed version is 
needed. First or second generation quarks are labelled q and top quarks are labelled t. 
Gluons are denoted as g and always represented by curly lines. Momenta and helicities 
are given in parentheses behind the label. For example, t{p, A) denotes a top quark with 
(four-)momentum p and helicity A. The generic scalar, vector and fermionic particles in 
each topology are given uppercase labels S, V, F, etc. For each topology we provide a 
list of MSSM particles that have to be substituted for the generic ones. For the various 
MSSM particles we use the notations from [8]. However, unless stated otherwise the 
generation indices /, J, etc. only run over the first two generations. The third generation 
quarks are written explicitly as t (= u^) and b (= d^). 

3.1 Tree-Level Diagrams 

The following prototype diagrams contribute to ti production at tree-level: 




Dqqbar_sV_tree 

V = 7,Z,g 



Dgg_sG_tree 
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Dgg_tF_tree* 



3.2 Self- Energy Diagrams 

For our calculation we only need to consider self-energy corrections for fermions and 
gluons. The SUSY self-energy corrections to fermion propagators only come from scalar 
particles: 

S{p-q) 

> ^ > ' > (3.1) 

P F{q) P 

The s-channel gluon propagator gets SUSY corrections from fermion and scalar loops: 



F{q) 5(g) 




OTJOT^ v , A* '"O'o'o'cr^ VwjOT^ v , M V .o.os.osix's.o.os.o. / v (3.2) 
p P \^ ^ p P P 

F{q-p) S{q-p) 

The two scalar self-energy diagrams always appear as a pair with the same coefficient 
and the same scalar particle in the loop. By inserting these self-energy corrections in 
individual lines of the tree level diagrams we obtain the following self-energy diagrams: 
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qiPa, Aa) 




t{pi, Ai 



OWO"0'OWOOOO> 



+ 




q{PaAa) t{pi,Xi 



Aft) 



Dqqbar_sV_xseSq 

(g,F,^) = (w^X°,t/.), (c^',X°,A), (w^X„A), (c^^x5,A), (w^A,A), (ci^A,A) 



^(^•a, Aa) 




A 



■0"0"0"0"0"0"0"0"0"0"0"0"0" 



q{pb,Xb) t{p2,X 




+ 




q{Pa,Xa) t{pi,\ 



■0"0"0"0"0'0"0"0"0"0"0"0"0"0" 



Afe) 




i(P2, A2) 



Dqqbar_sV_xseSt 



(F,5) = (x°,t/i), (x„A), (t,A?), (A,t/.) 




qiPa,Xa) t{pi,Xi) 




q{pb,Xb) t{p2,X2) 




Dqqbar_sG_iseF, Dqqbar_sG_iseS 
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Dgg_sG_xseSt 

(F,5) = (x°,?7,), (x„A), it,H^), {b,H^), {A,U,) 




Dgg_sG_iseF, gg_sG_iseS 

F = A , S = Ui,Di 




{F,S) = {x',,Ui), {t,Hf), (t,A?), {b,H^), {A,U,) 





Dgg_tF_xseSt* 

F = A , S = U,,D, 
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{F,S) = {xlU,), (x„A), {t,Hf), {b,H^), (A,t/,) 

In each diagram the hatched blob stands for one of the self-energy corrections from 
dSlD or (D- 
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3.3 Vertex Corrections 

The prototype vertex corrections for the qq tt amphtude are: 




Dqqbar_sG_vertSSq Dqqbar_sG_vertSSt 

{q,F,S)^ q^u',d' 

{u^,A,U,), {u',x',,U.), {u',Xj,D,), {S,F)^{A,Ui), {x'j,Ui), (x,-,A) 
(ci^A,A), {d',x'j,Di), {d',Xj,Ui) 
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For the gg — > tt amplitude we distinguish vertex corrections for s and t-channel 
diagrams. The corrections to the s-channel diagrams are: 




QsmiQsmiasiSiasia. 



g{Pb, Aft) 




Dgg_sG_vertSt 
(F,5) = (A,C/,), (i,^S), {b,H,) 




nFiq) 



Dgg_sG_vertSSt 
(5, F) = A), (C/,,x°), (A,X,) 




OOWO"0"0"0"OW 



*(P2, A: 




5(Po, Aa) 




0"0"0"0"0"0"0"0"0"0"ff 



*(P2, A 



5(P6, Ab) 




Dgg_sG_vertFg 

F = A 





ff(Pa, Aa) 



*(pi,Ai) 




9{Pb,h) 




Dgg_sG_vertSg 
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Dgg_sS_vertFg 
F = t,b , S = Hf,Al 




Dgg_sS_vertSg 

S,^Ui,Di , S-2^HlA\ 




Dgg_sS_vertSSg 
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The i-channel vertex corrections are: 
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3.4 Box Diagrams 



Si 



S2iq) 




Dqqbar_boxSS 

(g,5i,52,Fi,F2) = (w^t/„^7„A,A), {u^,U,,U,,A,xl), {u' ,U,,U,,xl 
(d^ A, Dj, A, Xfe), {d^: Ui, Uj, Xk: A) 



I 

I 



{A,B) 



I 
I 



q{pb,h) t{p2,X2) 

Dqqbar_fboxSS 

{d',A,A,Di,U,), {u',xl,A,Ui,Uj), 
{u',A,xl,Ui,U,), {d^,xl,A,D,,Uj), 
{d^,A,xlD„Uj) 



Fl 



52"^ N Si 



/ {A,B) 



^ ^2(9) 

^q{pb,h) t{p2,X2) 

Dqqbar_fboxSSx** 

(g, Fl, F2, Si, 52) = (m^ A, A, 

(d^A,A,A,[/,), («^x^,A,f/i,f/i), 
(^^^A,x^,C/^,C/,■), {d^,xlKDi,Uj), 
id',A,xlD„Uj) 



-^(PaiXa) 



I 

k 

I 



Dgg_boxFS* 

(F,5) = (A,t/,), 



9{Pa,K) 




I 

A 

I 





Dgg_boxSF* 

{S,F) = iU„A), {Ui,x'j), (D^Xj) 
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Dgg_boxSFx* 

(5,F) = (f/i,A) 



4 Cross Sections and Asymmetries 

In this section we present our results for the SUSY corrections to polarised tt production 
cross sections, which we calculated for each of the 10 Snowmass benchmarks detailed in 
[S]. To calculate the masses of the supersymmetric particles and run the couplings to 
the TeV scale we used the program SOFTSUSY by B. C. AUanach [6J. The decay widths 
of the MSSM Higgs particles were calculated with the program HDECAY by Djouadi, 
Kalinowski and Spira [7]. The Feynman rules for the MSSM vertices were taken from J. 
Rosiek's paper [8]. We compare our parton level cross sections with the results obtained 
in the leading log approximation [9]. Then we discuss our results for the total pp tt 
cross section and the double helicity asymmetries introduced in |10j . 



Let (jj denote the total cross section for the process i —>■ ti, where the initial state i 
can be a gluon pair (gg), a light up- type quark-antiquark pair (uu) or a light down type 
quark-antiquark pair (dd). We regard o"j as a function of the variable s = M^|, where 
Mtt is the invariant mass of the top-antitop pair. For each of these cross sections we 
have calculated the leading order contribution af'^ and the SUSY corrections due 
to the diagrams listed in the previous section. The SUSY corrections can be split into 
super-QCD (SQCD) corrections and super-electroweak (SEW) corrections. The SQCD 
corrections are of order 0{al) and the SEW correction of order O(aa^). Consequently 
the SEW corrections are one order of magnitude smaller than the SQCD corrections. 
We also define the ratios 

Figures O [6| and show a comparison of our "exact" ratios with the results obtained 
in the leading log approximation by Beccaria, Renard and Verzegnassi We have 

used Snowmass benchmark 5 to compute the exact cross section, but the observations 
stated here are true for any of the 10 Snowmass benchmarks. The only SUSY inputs 
in the leading log approximation are tan/5 and a universal SUSY mass scale Msusy- 
Sensible values for this scale lie anywhere between the mass of the lightest and the mass 
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Figure 5: SUSY corrections to the gg — > it cross section in tlie "exact" calculation and the 
leading log approximation for benchmark 5 of [5] . The width of the leading log graph reflects 
the uncertaincy due to the choice of the universal SUSY scale. 
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Figure 6: SUSY corrections to the uu — > ti cross section in the "exact" calculation and the 
leading log approximation for benchmark 5 of 0. The results for first and second generation 
up-type quarks are identical. The width of the leading log graph reflects the uncertainty due 
to the choice of the universal SUSY scale. 
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Figure 7: SUSY corrections to the dd — > tt cross section in the "exact" calculation and the 
leading log approximation for benchmark 5 of [5]. The results for first and second generation 
down-type quarks are identical. The width of the leading log graph reflects the uncertainty 
due to the choice of the universal SUSY scale. 



of the heaviest SUSY particle. The widths of the leading log graphs in figures El El and 
reflect this uncertainly. We see that, in the leading log approximation, the ratios ff^^"^ 
are proportional to log(s/M|ugY)- ^ ~ 10^ GeV the exact ratio runs linear with the 
same slope, but with a constant offset to the leading log graph. For very large centre 
of mass energies this offset becomes negligible. Therefore our results agree with the 
leading log approximation in the high energy limit. However, for collision energies that 
are achievable at the LHC we see that the leading log approximation fails to reproduce 
the exact results of the full SUSY calculation, even with the offset included. 

To obtain the pp — > ti cross sections, the parton level cross sections (Tj were folded 
with the CTEQ6L1 set of the CTEQ v6.51 parton distribution functions [IT]. For the 
proton-proton collision we assumed a centre of mass energy of 14TeV. Let da\^\^/ dMti 
denote the invariant mass differential cross section for producing a top quark with 
helicity Ai and an anti-top quark with helicity A2. Then we define 



dat, 



,ot 



dMu 
daiL 
dMu 
da pv 
dMn 



dMtt 

dMtt 
da^- 



+ 



+ 



da 


da+_ 
dMtt 
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For each combination we indicate the leading order and SUSY contributions by su- 



perscripts 'LO' and 'SUSY', respectively. The parity even combinations da, 
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Figure 8: Leading order results for the invariant mass differential cross section. 



daf^^^ are dominated by the SQCD corrections. However, for the parity odd combi- 
nation dap^^ the SQCD corrections are zero, since parity is conserved in super-QCD. 
For the asymmetries and the SUSY corrections we define the ratios 

dal2/Pv/dMu , , _ da^ot/lL/Pv/dMu 
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Figure 9: Leading order results for the invariant mass LL asymmetry. 
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Figure 10: SUSY corrections to the invariant mass differential cross section for the Snowmass 
benchmarks. The numbers in the legend refer to the labelling of the benchmarks in [5]. 

Figures El and [P| show the results for da^^/dMti and rYl{Mtt), respectively. Since 
there is no parity violation at leading order the ratio rpy is identically zero. 

Figures [TU and [HI show the ratios rf^^^^{Mtt) and rfl^^{Mtt) for each of the 10 
Snowmass benchmarks. We see that the SUSY corrections to the tt cross section 
can be as large as 10% of the leading order cross section, but typically only reach the 
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Figure 11: SUSY corrections to the invariant mass LL asymmetry for the Snowmass bench- 
marks. The numbers in the legend refer to the labelling of the benchmarks in [5]. 
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Figure 12: SUSY corrections to the invariant mass PV asymmetry for the Snowmass 
benchmarks. The numbers in the legend refer to the labehing of the benchmarks in [5]. 

5% level. In both plots we see "resonance peaks" located at the masses of the heavy 
and the pseudo-scalar Higgs {H^ and in the notation of [8J). They come from the 
scalar s-channel propagators in the diagrams labelled Dgg_sS_vertFg, Dgg_sS_vertSg 
and Dgg_sS_vertSSg in section [3 Note that these particles have negligible coupling 
to the light quarks extracted from the incoming hadrons, so that they arise from the 
gluon fusion contribution, in which the exchanged scalars are connected to the incoming 
gluons via a gluino or squark triangle. Moreover, due to the fermion triangle the sign 
of the contribution from diagram Dgg_sS_vertFg is opposite that of Dgg_sS_vertSg 
and Dgg_sS_vertSSg. This explains why we get downward-pointing resonance peaks 
for some of the benchmarks. Also note that, for all 10 benchmarks, the difference of 
the heavy and the pseudo-scalar Higgs masses is much smaller than their decay widths. 
Consequently we can only see two distinct peaks when these peaks have opposite signs. 
The kinks occurring between 1 and 2 TeV coincide, for each benchmark, with twice 
the gluino mass and can therefore be understood as a threshold effect due to the box 
diagrams Dqqbar_f boxSS and Dgg_boxFS. 

Figure [22 shows the SUSY corrections to the parity violating asymmetry for each 
of the 10 benchmarks. Here the resonance peaks are absent, because the diagrams 
Dgg_sS_vertFg, Dgg_sS_vertSg and Dgg_sS_vertSSg are parity-conserving. Further- 
more, the SUSY corrections to the parity violating asymmetry are one order of mag- 
nitude smaller than the corrections to the parity-even observables, because it only gets 
contributions from interferences of oder (^(aa^). 

By integrating the differential cross sections (14. 2 p over we obtain cross sections 
for producing tt pairs with arbitrary invariant mass. We define the cross sections dtot, 
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"tot 


— -O U O 1 
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"ll 


— D U o 1 


— -D U O 1 


la 
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+8.4 


91.99 


+1.75 


+0.063 


lb 


288.3 


+6.9 


79.50 


+1.72 


+0.020 


o 
Z 


zoo. 3 


—1.0 


72.84 


—0.34 


+0.125 


3 


289.5 


+6.9 


80.08 


+1.73 


+0.043 


4 


307.6 


+7.1 


85.62 


+1.64 


-0.018 


5 


332.1 


+19.8 


93.48 


+5.16 


+0.111 


6 


310.8 


+7.4 


85.74 


+1.68 


+0.053 


7 


284.9 


+3.5 


78.04 


+0.78 


+0.029 


8 


271.2 


-1.9 


74.08 


-0.64 


+0.035 


9 


263.6 


+6.9 


72.32 


+ 1.82 


+0.019 



Table 1: Numerical results for the integrated ti cross section and asymmetries. The numbers 
in the left column refer to the labelling of the Snowmass benchmarks in |5j. The superscripts 
'LO' and 'SUSY' indicate leading order results and SUSY corrections, respectively. The cross 
sections are given in pico-barns (pb). 



aiL and apv by 



I -J tot ILL I PV I . 

OtotiLLiPY = / dMu — , (4.4) 

J2mt U'J-VJ-tt 

where rrit is the top mass and Mpp = 14 TeV is the invariant mass of the proton-proton 
system. Again, the leading order and SUSY contributions are indicated by superscripts 
'LO' and 'SUSY', respectively. Table [I] summarises our results for these cross sections. 
For both, atot and ctll we see that the SUSY corrections typically make up 2% of the 
leading order results. However, they can be as big as 5% in the case of benchmark 5 
and as small as 0.3% in the case of benchmark 2. 

Experimentally it is often more convenient to parametrise the ti production cross 
section by the transverse momentum px of the top quark. For the transverse momen- 
tum differential cross section da\^\^/dpT we define the total differential cross section 
datot/dpx, the asymmetries daii/dpT and dapv/dpx and the ratios rYKpr), ''^YviPT)^ 
^tot'^i.PT): ^Tl^^^Pt) and r^^y^^px) in analogy to (14. 2 p and (14.31) . Our results for these 
quantities are shown in figures [I3| to [13 We note here that the "resonance peaks" 
and "troughs" from the thresholds for scalar particle exchange are smoothed out by 
the phase-space integration which means that Mti is a far better variable in which to 
analyse the data in order to extract information on the SUSY parameter set, although 
we note that some of the benchmarks give rise to an enhancement of the differential 
cross-section of up to 7% at large pr- 
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Figure 13: Leading order results for the transverse momentum differential cross section. 
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Figure 14: Leading order results for the transverse momentum LL asymmetry. 
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Figure 15: SUSY corrections to the transverse momentum differential cross section for tlie 
Snowmass benchmarks. The numbers in the legend refer to the labelling of the benchmarks 
in ^. 
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Figure 16: SUSY corrections to the transverse momentum LL asymmetry for the Snowmass 
benchmarks. The numbers in the legend refer to the labelling of the benchmarks in [5]. 



28 



0.003 




-0.001 I ' ^ ' ^ ' 1 ' ^ ' 1 

100 200 300 400 500 

Pt [GeV] 



Figure 17: SUSY corrections to the transverse momentum PV asymmetry for the Snowmass 
benchmarks. The numbers in the legend refer to the labelhng of the benchmarks in [5]. 

5 Conclusions 



We have calculated the complete MSSM corrections to the cross-section for t — i produc- 
tion at LHC The calculation has been set up in terms of a set of prototype Feynman 
graphs for the partonic-level helicity matrix-elements. In this way the prototypes can 
be used for any set of SUSY parameters by inserting the corresponding couplings and 
internal masses into the prototypes. 

We have obtained numerical results for the ten Snowmass benchmark sets using the 
CTEQ PDF's. We find a considerable variation of the effects of the one-loop SUSY 
corrections between the various benchmarks. The benchmark giving the largest correc- 
tion is benchmark 5, which is a super-gravity model with small tan/? = 5 and a large 
negative tri-linear coupling, Aq = —1000 GeV. These large corrections can be under- 
stood from the fact that this large tri-linear coupling generates a light stop mass (258 
GeV) thereby enhancing graphs involving a stop mass inside the loop. This gives an 
enhancement of 6% in the total production cross- sect ion. 

Whereas the corrections for the other benchmarks are somewhat smaller, they are 
usually around 3% and therefroe comaparable to the weak corrections calculated by 
Bernreuther et. al. [12] and Kuhn et. al. [13]. Note that whereas the weak corrections 
reported in [12], [13] decrease the prediction for the cross-sections, the SUSY corrections 
have a positive sign for most of the benchmarks considered. 

Statistically, we expect these events to be easily detectable given the anticipated yield 
of order 10^ events over the period of running of LHC. We have found similar corrections 
in the asymmetry ratios defined in in Eqs. (14.21) and (14. 3p . For these asymmetries we also 
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expect cancellation of systematic errors arising from uncertainties in incoming parton 
fluxes, so that these corrections of order 3% would exceed the statisical errors by a factor 
of~ 100. 

Given corrections of such significance, it is reasonable to assume that corrections 
in the differential cross-sections will also be detectable (provided sufficiently large bins 
are taken). We have therefore plotted the differential cross-sections with respect to the 
invariant mass Mft of the t — i system and also with respect to the transverse momentum 
Pt of th t-quark. In the former case, the differential cross-sections display an interesting 
structure with peaks and/or troughs corresponding to thresholds for scalar particle 
exchanges in the gluon fusion process. 

We have also determined the SUSY contribution to the parity odd hehcity asym- 
metry. This receives only contributions from the supersymmetric partners in the weak- 
interaction sector, which are suppressed relative to the SQCD corrections by 0{aw / (^s)- 
It would appear, therefore that even for benchmark 5, which produces the largest cor- 
rections, such parity violating asymmetries will be too small to observe. 
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